Controlled generation and steering of spatial gap solitons 
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We demonstrate the first fully controlled generation of immobile and slow spatial gap solitons in 
nonlinear periodic systems with band-gap spectra, and reveal the key features of gap solitons which 
distinguish them from conventional counterparts, including a dynamical transformation of gap soli- 
tons due to nonlinear inter-band coupling. We also predict theoretically and confirm experimentally 
the effect of anomalous steering of gap solitons in optically-induced photonic lattices. 

PACS numbers: 42.65.Tg,42.65.Jx,42.70.Qs 



Nonlinear wave self-action plays an important role in 
many physical systems ranging from water waves to opti- 
cal beams and Bose-Einstein condensates (BECs). Non- 
linearity can suppress wave spreading leading to the for- 
mation of solitary waves which, due to their robustness, 
can find many applications, e.g. in long-distance com- 
munication systems. Almost two decades ago, it was 
suggested that the systems with periodically modulated 
parameters can support a novel type of solitons-^ap soli- 
tons^ which may exist in band gaps of the linear spec- 
tra due to strong wave scattering and coupling between 
the forward and backward propagating waves 0, • It 
was predicted that gap solitons may exist in different 
nonlinear periodic structures including fiber Bragg grat- 
ings , two- and three-dimensional photonic crystals ^ , 
and BEC in optical lattices 0- 

According to the theoretical predictions, gap solitons 
have many unique properties, which distinguish them 
from conventional solitons (i) they can form in 

both self-focusing and self-defocusing nonlinear media; 
(ii) soliton velocity and dispersion are nontrivially mod- 
ified; (iii) the gap solitons become unstable and undergo 
dynamical transformations above a critical amplitude 
due to the inter-band resonances ^ ^7J . For gap solitons 
observed experimentally in fiber Bragg gratings |^ and 
waveguide arrays 9], only some reduction of the soliton 
velocity was observed. This has prevented an experi- 
mental study of dispersion properties of gap solitons and 
associated band-gap effects, since such effects become ev- 
ident for slow and immobile solitons only. 

In this Letter, we demonstrate the first fully controlled 
generation of spatial gap solitons where the soliton power 
and velocity may be independently varied. In particular, 
we observe experimentally immobile gap solitons. We 
demonstrate the gap-soliton steering, reveal their unusual 
mobility, and observe, for the first time to our knowledge, 
the gap-soliton dynamics due to instability development. 
Our experimental results are in excellent agreement with 
the key theoretical predictions, underlying generic prop- 



erties of gap solitons with implications to a variety of 
nonlinear periodic systems beyond the field of optics. 

We investigate the gap-soliton formation in photonic 
lattices created by interfering ordinary polarized laser 
beams in a biased photorefractive crystal. Discrete soli- 
tons in such lattices were observed recently 0, ^ ; this 
allows us to compare two types of solitons emphasizing 
the intriguing properties of gap solitons. 

First, we determine the optimal conditions for gap- 
soliton generation, using the normalized paraxial equa- 
tion for the probe beam electric field envelope E{x^z), 
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■r{x,\E\'')E = 0, 



(1) 



where x and z are the transverse and propagation coordi- 
nates, normalized to the characteristic values xq and zq, 
respectively, D = zoA/(47rnoXo) is the beam diffraction 
coefficient, no is the average medium refractive index, 
and A is the vacuum wavelength. The optically induced 
change of the refractive index is 0, 



J^ix,\Ef) = -jih+Ipix) + \Ef)- 



(2) 



where is the constant dark irradiance, Ip{x) is the 
two-beam interference pattern which induces the periodic 
lattice with a period 



Ip{x) 



Ig cos^ {irx I d) , 



(3) 



and 7 is a nonlinear coefficient proportional to the ap- 
plied DC field. In order to match our experimental condi- 
tions, we use below the following parameters for the the- 
ory and numerical calculations: A = 0.532 /im, no = 2.4, 
Xq = 1 /im, zq = 1mm, d = 22.2, /5 = 1, = 1, 
7 = 5.31, the crystal length is L = 15 mm. We note 
that the model Eq. Q is very general and appears in 
other areas of physics. It can describe, in particular, the 
BEC dynamics in optical lattices where the matter-wave 
gap solitons are expected to have similar properties to 
optical counterparts. 
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FIG. 1: (a) Dispersion of Bloch waves in an optically-induced 
lattice; the spectrum bands are shaded, (b) Bloch-wave pro- 
files (solid) and leading-order Fourier components (dashed) 
superimposed on top of the normalized refracted index pro- 
file of the periodic lattice (shown with shading) for different 
gap edges indicated by the arrows. 



Existence of spatial bright solitons is closely linked to 
the structure of the linear wave spectrum. In periodic 
lattices, the spectrum is composed of bands correspond- 
ing to the propagating Floquet-Bloch modes^ which are 
separated by gaps where the wave propagation is forbid- 
den. The Floquet-Bloch modes are solutions of linearized 
Eq. (QJ in the form E^^n{x,z) = i/j^^ni^) ex.p{ihzx/d -\- 
^/^/^,n^), where P^^n and tz are the Bloch-wave propa- 
gation constant and wavenumber, respectively, and the 
index n = 1,2, .. . marks the order of the transmission 
band. The Bloch wave '0,^,n(^) has the periodicity of the 
photonic structure, tpK,n{x -\- d) = i/j^^ni^)^ and this con- 
dition defines the dispersion relation 0^,71^ as sketched in 
Fig. da). The open regions in this plot mark band gaps 
where Im^^ 7^ 0: the top one exists due to total internal 
reflection and extends to /3 +00, whereas lower gaps 
have a finite width and appear due to Bragg scattering 
from the periodic structure. 

In order to underline the essential physics of the soliton 
formation in periodic lattices, first we consider the small- 
amplitude limit. We seek solutions of Eq. (QJ near the 
band edge in the form of modulated Bloch waves, E = 
(f{x^ z)tlj^^ri e^'p{i/3^^nZ-\-ii^x/d)^ and derive the nonlinear 
Schrodinger equations for the slowly varying Bloch-wave 
envelope 12], 
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ir4 + r|^IV = o. 



(4) 



Here Vg = —ddl3/dn\^ is the group velocity, r] = 
—d^ i5 1 di<?\^ is the effective diffraction coefficient, 

and F = l^{T{x,\^^^J')-T{x,^)\\^^^J'dx is the 




FIG. 2: Numerical results for a gap soliton family: Power 
(top) and width (middle) vs. the propagation constant. Bot- 
tom: soliton profiles (solid) corresponding to the marked 
points (a) and (b) in the upper plots; shadings mark the 
grating minima. Arrows illustrate schematically the direc- 
tion of the input beams, which interference pattern (shown 
with dashed line) approximates the soliton profile. 



effective nonlinear coefficient, where we assume the 
normalization ^'^\\\)t^^n\^dx = 1. Equation pos- 
sesses localized solutions for bright solitons, ip = 
Y^2p/Fsech {x — Vgz)^/p/r] exp(zpz), provided r^F > 
and p is a free parameter. Thus, for self- focusing medium 
nonlinearity (F > 0), the gap solitons appear near the 
lower-gap edges of the bandgap spectrum. 

Next, we find soliton solutions of the full model 
Eqs. (QJ-©, shown in Fig. [2j and identify a number 
of important differences between the conventional and 
gap solitons which appear when the propagation con- 
stant is shifted deep inside the gap, where the envelope 
approximation is no longer valid. Specifically, the 
power of gap solitons is bounded from above [Fig.[2ftop)], 
and their width is bounded from below [Fig. [21 middle)]. 
These limitations appear because the Bragg-reflection 
gap always has a finite width, in contrast to the semi- 
infinite gap where conventional (or discrete) solitons are 
localized. 

Controlled experimental excitation of spatial gap soli- 
tons can be realized if the modulated Bloch-wave profile 
is properly matched at the input. Since the Bloch waves 
are periodic, they can be decomposed into the Fourier 
series, tlj,<c,n — ^n{f^ + 27rm) e'x.-p[ix{hi + 27Tm)/d). In 
Fig. H^b) we show the characteristic profiles of the Bloch 
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waves, and also plot the contribution from the leading- 
order Fourier components (dashed lines). We find that 
in the leading order ?/^o,i = Ci(0) + . . ., therefore lattice 
solitons in the semi-infinite total internal refiection gap 
can be generated by a single incident beam, as was real- 
ized in earlier experiments for arrays of weakly coupled 
optical waveguides. On the contrary, the Bloch waves in 
the Bragg-reflection gap are composed of the counter- 
propagating waves, e.g. ip^^n = Cn(7r) exp(z7rx/(i) + 
Cn{—7T) exp(— ZTTx/d) + . . . for n = 1, 2. Therefore, spa- 
tial gap solitons can be generated by using two Gaussian 
beams which are tuned to the Bragg resonance and have 
opposite inclination angles, as was originally suggested 
in Ref. 13] and developed further in Ref. 7]. Therefore, 
we consider the input electric field in the form. 



Eo{x) 



-(x-Xe)V^%os[7r(x-x,)/d], (5) 



where the exponential term approximates the gap soliton 
envelope, w being the width of the input beams, and the 
interference term approximates the Bloch- wave profile, 
with the shift Xs depending on the relative phase differ- 
ence between the two beams, and we choose = d/2 
and w = 55. When the interference maxima are at the 
minima of the refractive index profile, the Bloch mode is 
excited at the lower edge of the Bragg-refiection gap, and 
the input pattern can match very closely the gap-soliton 
profile, as shown in examples of Figs. [2fa,b). 

We investigated numerically the dynamics of the two- 
beam mutual focusing and the gap-soliton generation by 
simulating the model Eq. with the input condition Q 
chosen to match the Bloch-wave profile at the lower gap 
edge. The beams diffract at a low input power [Fig.[3Ja)], 
whereas an immobile gap soliton forms when the input 
power is increased [Fig.IH^b)]. We note that the required 
power depends on the input beam width, as follows from 
Fig. [21 and the minimum soliton width defines a fun- 
damental limitation of the degree of two-beam mutual 
focusing. Indeed, as the power is increased, an instabil- 
ity rapidly develops through a resonant excitation of the 
first band, and subsequent formation of a quasi-periodic 
breather [Fig. Eljc)]. Similar breathing states were re- 
cently observed in waveguide arrays |l4j . 

Our experiments were performed in a 15mm- long 
Strontium Barium Niobate (SBN:60) crystal externally 
biased along the crystalline c-axis. The experimental 
setup is similar to that discussed in Ref. [Il| with the 
difference that the extraordinary polarized probe beam 
was split into two parts which were later focused by cylin- 
drical lens and made to overlap at the input face of the 
crystal. The angle between these two beams was set to 
twice the Bragg angle, such that the periodicity of the in- 
terference pattern is equal to that of the lattice (22 /im). 
In our case, such value of the periodicity allows for a rel- 
atively wide gap in the transmission spectrum, as shown 
in Fig.nja). The relative phase between the probe beams 
was tuned such that a symmetric interference pattern is 
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FIG. 3: Numerical results. Dynamics of the Bloch waves 
excited through two-beam interference: (a) linear diffraction 
at low power (/o ^ 0), (b) excitation of a gap soliton in the 
nonlinear regime (/o = 0.048), (c) beam breakup and the 
formation of a quasi-periodic breathing state at higher powers 
(/o = 0.28). Left: variation of intensity along the propagation 
direction; Right: beam profiles at the crystal output {z = 
15mm) normalized to Iq. 



obtained, as shown in Fig. ^top). The relative position 
between this pattern and the lattice could also be con- 
trolled, by changing the relative phase between the two 
lattice forming beams !TT|. The input width of the over- 
lapping probe beams is w = 55/im (65/im FWHM), while 
they are fully separated at the output for zero bias field. 
When electric field of 5000 V/cm is applied to the crys- 
tal, the index grating (optical lattice) is formed and the 
probe beams excite Bloch waves in the first or the second 
band, depending on the relative lattice position. 

The experimental results shown in Fig. ^represent re- 
markably all theoretical predictions. First, we align the 
interference maxima of the input beams with the min- 
ima of the induced lattice at the input face of the crys- 
tal and record the beam profiles at back face for sev- 
eral input powers [see Fig. ^Jjleft column)]. The out- 
put intensity is exactly zero at the maxima of the in- 
dex grating. The intensity maxima are out-of-phase, as 
confirmed with inter ferometric measurements, and posses 
a double peak structure located at the minima of the 
grating. At low powers, the output intensity pattern is 
broad and corresponds to the Bloch waves at the lower 
gap edge. The beam is exactly centered between two 
unperturbed output beams [dashed curve in Fig. ^left 
column, top plot)] measured for zero voltage. At higher 
intensity {lo/Ig = 0.04), we observe mutual focusing, 
and self-trapping into a spatial gap soliton is observed 
at lo/Ig = 0.14 [see Fig. EJleft column)]. The gener- 
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FIG. 4: Experimental results for two-beam interaction show- 
ing the intensity profiles at the crystal input (top) and output 
for varying beam power: Left: mutual focusing and gap soli- 
ton formation when interference maxima are aligned with the 
lattice minima; Right: self-defocusing when interference max- 
ima of the input beams are at the lattice maxima. Dashed 
curves represent the beam profiles at the indicated beam in- 
tensity when the grating is erased. 
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FIG. 5: (a) Schematic demonstration of anomalous gap- 
soliton steering induced by a tilt of the lattice; dashed line 
shows the propagation direction of a gap soliton; solid lines 
indicate the directions of the input beams. (b,c) Output soli- 
ton profile for a lattice tilt in the direction of larger x by 20% 
of the Bragg angle with respect to normal; dashed lines show 
the beam profiles when the lattice is absent. 



ated gap soliton has zero transverse velocity, and is cen- 
tered between the two output beams which separate if 
the grating is erased (dashed curve). As predicted theo- 
retically, the effect of the mutual focusing is limited and 
at higher intensities {lo/Ig = 0.25) the beam disinte- 
grates [see Fig. 0Jleft column, bottom plot)] while its 
profile becomes asymmetric due to the diffusion contri- 
bution to the photorefractive nonlinearity. On the other 



hand, when we align the interference maxima of the in- 
put beams to the lattice maxima, the excited Bloch wave 
corresponds to the upper gap edge [see Fig. Q and they 
experience anomalous diffraction {r] < 0) leading to self- 
defocusing as the power is increased [Fig. fright)]. 

Finally, we study mobility of spatial gap solitons and a 
possibility to vary their transverse velocity. Experimen- 
tally, the soliton motion is induced by tilting the lattice 
by 20% of the Bragg angle thus introducing a lateral shift 
of the induced waveguides by 16/im at the output [the 
lattice is shifted to the right in Fig.ISfa)]. Results of the 
experimental observations and the numerical simulations 
are presented side-by-side in Fig. IStb,c) show that the 
generated gap solitons move to the left when the grat- 
ing is tiled to the right. In experiment, on the right-hand 
side of the gap soliton, a small contribution from the first 
band is observed, which appears due to asymmetry in the 
initial excitation and inhomogeneities of the lattice. This 
anomalous steering behavior occurs because the spatial 
group- velocity dispersion for gap solitons is almost three 
times larger compared to a homogeneous crystal. This is 
in sharp contrast with conventional lattice (or discrete) 
solitons, that tend to propagate along the lattice 15]. 

In conclusion, we have demonstrated experimentally 
the first fully controlled generation of spatial gap soli- 
tons in optically-induced periodic photonic lattices and 
observed novel effects such as anomalous steering of gap 
solitons and the limitation of the two-beam mutual fo- 
cusing through inter-band coupling. We believe our re- 
sults can be useful for the study of nonlinear effects in 
photonic crystals and nonlinear dynamics of the Bose- 
Einstein condensates in optical lattices. 
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